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Abstract 

We studied the acceleration mechanisms of chromospheric jets associated with 
emerging flux using a two dimensional magnetohydrodynamic (MHD) simulation. 
We found that slow mode shock waves generated by magnetic reconnection in the 
chromosphere and the photosphere play key roles in the acceleration mechanisms of 
chromospheric jets. An important parameter is the height of magnetic reconnection. 
When magnetic reconnection takes place near the photosphere, the reconnection out- 
flow collides with the region where the plasma beta is much larger than unity. Then 
the plasma moves along a magnetic field. This motion generates a slow mode wave. 
The slow mode wave develops to a strong slow shock as it propagates upward. When 
the slow shock crosses the transition region, the transition region is lifted up. As a 
result, we obtain a chromospheric jet as the lifted transition region. When magnetic 
reconnection takes place in the upper chromosphere, the chromospheric plasma is 
accelerated due to the combination of the Lorentz force and the whip-like motion 
of magnetic field. We found that the chromospheric plasma is further accelerated 
through the interaction between the transition region (steep density gradient) and a 
slow shock emanating from the reconnection point. This is an MHD effect which has 
not been discussed before. 

Key words: Sun: chromosphere — Sun: photosphere — Magnetic reconnec- 
tion — Shock waves 



1. Introduction 

Collimated ejections of chromospheric plasma into coronal heights have been studied for 
many years. Spicules are jet-like structures seen in the limb of the Sun (Sterling, 2000; Anan 
et al., 2010). Surges are chromospheric jets seen dark in Ha images and their heights reach up 
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to ~ 200 Mm (Roy, 1973; Yoshimura et al., 2003). Recent observations have found the tiny 
chromospheric jets with the loop structures at their foot-points (Shibata et al., 2007). They 
are called the chromospheric anemone jets. Their typical height is 1 — 4 Mm (Nishizuka et al., 
2011). Figure 1 shows an example of the jets observed by the Solar Optical Telescope (Tsuneta 
et al., 2008) on board the Hinode satellite (Kosugi et al., 2007). 

Chromospheric jets are often observed in emerging flux regions and around sunspots. 
They are cool (~ 10 4 K) compared to coronal plasma (> 10 6 K). Many chromospheric jets 
disappear in the images of chromospheric lines before returning to the chromosphere. This may 
be due to some heating process or rapid expansion of the chromospheric material of jets (the 
rapid expansion is shown in this paper). Therefore for correct interpretations of observations 
we should understand how chromospheric plasma is ejected. 

Trigger of chromospheric jets can lie at various height in the solar atmosphere. X-ray 
observations found coronal hot jets (X-ray jets) associated with magnetic loops (Shibata et al., 
1994), which implies that the energy release process takes place in the corona. Multi- wavelength 
coronal observations have shown that some cool jets (surges) coexist with hot jets (Schmieder et 
al., 1995; Canfield et al., 1996; Chae et al., 1999), which implies that the trigger is located in the 
lower corona or in the upper chromosphere. It is known that some surges are associated with 
the Ellerman bombs or moustaches (Ellerman, 1917; Rust, 1968), where the Ellerman bombs 
are sudden energy releases in the low chromosphere or in the photosphere and their emissions 
are bright at both wings of the Ha line (Kitai, 1983; Fang et al., 2006; Watanabe et al., 2008; 
Matsumoto et al., 2008). This implies that they are triggered in the low chromosphere or in 
the photosphere. 

Using MHD simulations, Yokoyama & Shibata (1996) studied the acceleration mecha- 
nism of surges. They simulated the evolution of an emerging flux in a uniform coronal magnetic 
field. They obtained both hot and cool jets (which can correspond to the X-ray jets and surges, 
respectively) simultaneously as results of magnetic reconnection between the emerging flux and 
the pre-existing coronal field, where magnetic reconnection is a physical process in which a mag- 
netic field in a highly conducting plasma changes its topology due to finite resistivity. Their 
results are consistent with multi- wavelength observations. They pointed out that the hot jets 
and the cool jets are accelerated mainly due to the gas pressure gradient of the heated plasma 
and the Lorentz force, respectively. Note that the chromospheric evaporation due to the heat 
conduction is also important for the acceleration of the X-ray jets (Shimojo & Shibata, 2000). 
More elaborated modeling for coronal jets has been performed by other authors (Isobe et al., 
2006; Pariat et al, 2009; Archontis et al, 2010; Jiang et al, 2012). Martmez-Sykora et al. (2011) 
proposed that strong Lorentz force due to large magnetic field gradient squeezes chromospheric 
plasma and produces spicule-like jets moving at a local sound speed. Moreno-Insertis et al. 
(2008) compared a three-dimensional numerical experiment with an observed coronal jet and 
showed that their numerical results are consistent with the observations. 
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Fig. 1. (a) Call H broadband filter image of the active region NOAA 10935. (b) Enlarged image of the 
chromospheric anemone jet indicated by the white arrow in (a), (c) Left: Time-sequenced images of the 
chromospheric jet (Call) shown in (b). Right: Running difference images of the jet. Cross marks indicate 
the heads of the apparent horizontal motion along the loop structure and plus marks indicate the heads 
of the vertical jet. 1 arcsec~730 km. 



When magnetic loops expand in the solar atmosphere, magnetic reconnection can take 
place at various heights. Now we should note that in the chromosphere the gravitational 
stratification is much stronger than in the corona. As a result, the plasma beta varies from 
much larger than unity to less than unity. Therefore, the physical processes associated with 
magnetic reconnection can be different at various heights in the solar atmosphere. 

Can the lower atmospheric reconnection produce tall (~1-10 Mm) jets? For example, 
observations have suggested that chromospheric anemone jets are considered as results of mag- 
netic reconnection in the chromosphere (Morita et al., 2010; Singh et al., 2011). Their speed 
is comparable to the expected local Alfven speed 10 km s _1 ) in the upper chromosphere. 
Suppose that magnetic reconnection takes place in the upper chromosphere. Then it is highly 
possible that jets are reconnection outflows accelerated due to the whip-like motion of magnetic 
field, considering that the outflow speed is almost the same as the local Alfven speed. Next, 
suppose that magnetic reconnection occurred in the middle or low chromosphere where the 
plasma beta is expected to be larger than unity. Then it is difficult to regard observed jets as 
reconnection outflows because the Alfven velocity is lower than the observed velocity. 

Let us estimate the maximum height which reconnection outflows can reach. Suppose 
that the kinetic energy is converted into the potential energy. Then we obtain pV\j2 ~ pgHj et) 
where p, V^, g and Hj et are the density, the Alfven velocity, the gravitational acceleration and 
the height of the jet, respectively. We obtain Hj et <~ (p/ pg)(B 2 / (87179)) ~ H p /(3, where p, £>, 
H p and (3 are the gas pressure, the magnetic field strength, the pressure scale height and the 
plasma beta in the chromosphere, respectively. Suppose, for example, /3 ~ 1. Then the height of 
the jet will be H ~ 200 km, which is much shorter than the typical height of observed jets (e.g. 
1 — 4 x 10 3 km for the chromospheric anemone jets). Therefore we need to consider another 
scenario in which only a fraction of plasma is accelerated. 

A promising model is the acceleration mechanism through the interaction between a 
slow shock and the transition region (Osterbrock, 1961; Shibata et al., 1982, 2007; Suematsu et 
al., 1982; De Pontieu et al., 2004; Heggland et al., 2007). The transition region can be launched 
when a slow shock pass through it. The ejected chromospheric plasma behind the contact 
discontinuous layer is regarded as chromospheric jets. Shibata et al. (2007) pointed out that 
slow mode waves can be generated by magnetic reconnection between small emerging flux and 
the pre-existing ambient field. Then they proposed that the slow mode waves become shocks 
for some reason and produce chromospheric jets. Magnetic reconnection can generate waves of 
various kinds (Isobe et al., 2008; Nishizuka et al., 2008; He et al., 2009; Kigure et al., 2010). 
The wave can become shocks through non-linear processes. Carlsson & Stein (1997) compared 
observations with simulations and found that slow shocks are ubiquitously generated in the 
chromosphere. As supported in the above studies, the slow shock acceleration mechanism has 
the potential to systematically account for observational characteristics of chromospheric jets. 

The slow shock acceleration mechanism has been studied by many authors. Suematsu 
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et al. (1982) and Shibata et al. (1982) studied the mechanism using one-dimensional hydro- 
dynamic simulations. Their results are both qualitatively and quantitatively consistent with 
observational properties of spicules and Ha surges. They also proposed the possibility that the 
Ellerman bombs generate slow shocks and lead to surges. Saito et al. (2001) showed that the 
Alfven waves also produce slow shocks due to a non-linear coupling. Although both slow and 
fast mode waves are generated through a non-linear process of the Alfven waves, their results 
suggest that slow mode waves play more important roles in generation of jets than fast mode 
waves. Some authors discussed the slow shock mechanism by including other effects (Sterling et 
al., 1993; Heggland et al., 2007; Martmez-Sykora et al., 2009; Murawski et al., 2011). Hansteen 
et al. (2006) and De Pontieu et al. (2007) discussed the relation between chromospheric jets 
and the slow mode shocks generated by convective flows and global p-mode oscillations. 

In the previous studies, waves are generated in the chromosphere or in the photosphere, 
artificially (e.g. Shibata et al., 1982), by using observational data of the photospheric motions 
(e.g. Matsumoto et al., 2008), and convective flows in a self-consistent manner (e.g. Hansteen 
et al., 2006). Although Martmez-Sykora et al. (2009) suggested the possibility that sudden 
energy release due to reconnect ion in the low atmosphere can lead to chromospheric jets, they 
did not identify the connection. The relation between the acceleration of chromospheric jets 
and magnetic reconnect ion has been still unclear. 

Even if magnetic reconnection takes place in the upper chromosphere, the interaction 
between slow shocks and the transition region can play an important role in accelerating chromo- 
spheric jets. The standing slow shocks will emanate from a reconnection point if the resistivity 
is localized in some way (the Petschek-type reconnection; Petschek (1964)). For a detailed 
discussion on the possibility of the localization of the resistivity, see Section 4.5. Suppose that 
the Petschek-type reconnection occurred in the upper chromosphere. Then the reconnection 
inflow advects the transition region to the reconnection region, particularly to the slow shocks. 
As a result, the slow shock can cross the transition region. The interaction can change the flow 
structure of the reconnection outflow. Now we have less knowledge of the interaction processes 
in MHD shock cases. Therefore the interaction makes it difficult for us to correctly understand 
the dynamics of jets. 

To self-consistently account for the acceleration mechanisms of chromospheric jets associ- 
ated with emerging flux, we performed a two-dimensional (2D) MHD simulation. In particular, 
we focused on the relation between magnetic reconnection and chromospheric jets. Magnetic 
reconnection occurred between emerging magnetic loops and the pre-existing field in the chro- 
mosphere and at the foot-points of the emerging flux. Slow shocks are generated by magnetic 
reconnection in the lower atmosphere (near the photosphere) and upper chromosphere. The 
slow shocks accelerate the chromospheric plasma along a magnetic field, producing chromo- 
spheric jets. The key process in this paper is the interaction between slow shocks and the 
transition region. An important point we found is that the reconnection height determines the 
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acceleration mechanism of chromospheric jets. To investigate the interaction process in detail, 
we also performed numerous 1.5D simulations, as presented in Appendix. 

We note that evolution of emerging flux in a 3D space is similar to that in a 2D space. 
For emerging process, compare Shibata et al. (1989) (2D) and Matsumoto et al. (1993) (3D) 
and for Ellerman bombs, Isobe et al. (2007) (2D) and Archontis & Hood (2009) (3D). In 3D 
we will find a wider variety of physics than in 2D, including the interchange instability (Isobe 
et al., 2006) and magnetic reconnection between not perfectly anti-paralleled field lines (Jiang 
et al., 2011; Nakamura et al., 2012). Here we focus on fundamental physics in a 2D space. 

In Section 2, we introduce the numerical model adopted here. In Section 3, we present 
results of our 2D numerical simulation. Here we separately describe the acceleration processes 
of the two jets Jet A and B to clarify that the height of magnetic reconnection determines the 
acceleration mechanism. To investigate the acceleration process of Jet B in detail, we performed 
numerous 1.5D simulations, which is summarized in Appendix. In Section 4 we present a 
schematic diagram summarizing the acceleration mechanisms of chromospheric jets. Here we 
also give implications for observations and discuss the validity of the adopted assumptions. 

2. Numerical Setup 

2.1. Basic Equations and Assumptions 

We solve the two-dimensional, resistive and compressible MHD equations in Cartesian 
coordinates (x,z), where the z-direction is antiparallel to the gravitational acceleration. We 
assume that the medium is an inviscid perfect gas with a ratio of the specific heats of 7 = 5/3. 
We include the radiative cooling effect as described below but neglect the heat conduction (For 
the discussion on this validity, see section 4.6). The unit of length, velocity, time and density 
in the simulation are H , C s o, H /C s o = r and p 0) respectively, where H = k^T^/{mg^) is the 
pressure scale height (/cb is the Boltzmann constant and m is the mean molecular mass), C s o 
the sound speed and p the density at the base of the photosphere {z — 0). The gas pressure, 
temperature, magnetic field strength and energy are normalized by those units, i.e., p = PoC^, 
T = 771(7^0/(7^), ^0 = (PoCso) 1 ^ 2 an d E = PqCsqHq, respectively. The gravity is normalized 
by g = (7 2 / (^Hq). The values of the normalization units are given as the typical values in the 
photosphere: H = 170 km, C s o = 6.8 km s _1 , r = H /C s q = 25 s and p = lAx 10 -7 g cm -3 . 
Thus we obtain p = 6.3 x 10 4 dyn cm" 2 , T = 5,600 K, B = 250 G and E = 3.1 x 10 26 erg. 

The basic equations are as follows: 
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P=—pT (4) 

rn 

— -cVxE = (5) 

E = r]J--xB (6) 

c 

J = ^-\7xB (7) 

47T 

Here, g = (0, — go) is the gravitational acceleration, rj is the resistivity and J is the current 
density. T coo i ing represents the radiative cooling time given as a function of z. We include the 
radiative cooling effect in a simple manner: we assume that the radiative cooling time varies 
linearly in the photosphere and chromosphere and constant in the corona. The cooling term 
is only applied to the regions where the pressure fluctuation is positive compared to the initial 
pressure (therefore we do not apply the radiative cooling to the expanded regions) . We do not 
consider the radiative cooling in the convectively unstable zone {z < 0). 

f oc O<0) 

T C ooling{z) = < (j cl ~ T c0 )z/z tT + T c0 (0 < Z < Z tr ) (8) 

{ r cl [z > Z t r), 

where r c0 (= 0.04r =l s) and r c i(= 40r = 1,000 s) are the radiative cooling time at bottom and 
top of the chromosphere, respectively. Here we roughly mimic the average cooling time given 
in VAL-C model (Vernazza et al., 1981). z tr is the height of the transition region. We also 
assume the anomalous resistivity model in the simulation box. 

for Vd < v c 

a(v d /v c -l) 2 for v d >v c 



n = { , , M 2 r ( 9 ) 



Vd = \J\/p (10) 
Vc = Vci + (v C 2-v cl )^- tanh^^—^^+1 J, (11) 

where a(— 0.01Hq/t) : v d and v c are a constant, the (relative ion-electron) drift velocity and the 
threshold above which the anomalous resistivity sets in, respectively. v&(= 50C s o) and v c2 (= 
3,OOOC s o) are the threshold in the chromosphere and the corona, respectively and w c = 0.5H . 

2.2. Initial and Boundary Conditions 

We simply model the solar atmosphere: a super-adiabatically stratified layer representing 
the upper convection zone {Z^ n < z < 0), an isothermal cool (T = T ) layer representing the 
photosphere/chromosphere (0 < z < z tr ), and an isothermal hot (T = 150T ) layer representing 
the corona [z > z tr ). We take z = to be the base height of the photosphere. 
The initial distribution of temperature is given as 

dT 



T(z) = T pho - a 



dz 



z (Z min <z<0) (12) 
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Fig. 2. Schematic picture of the initial condition. Sold lines are the magnetic field lines. The computa- 
tional box is in the range — X m i n < x < X max and Z m [ n < z < Z max . 



for convectively unstable zone and 

T(z) = T pho + (T cor — T pho ) j- 



tanh 



(■- 



\ Wtt 



1 



} (*>0) 



(13) 



for the upper atmosphere, where T p h and T cor are the temperatures in the photo- 
sphere/chromosphere and the corona, respectively. We took z tr = 10H (1700km), w tr = 
0.2H (34km). \dT/dz\^ = (7 — l)/7 is the adiabatic temperature gradient and a is a dimen- 
sionless constant. The layer is convectively unstable in the case that a > 1. We used a = 1.5 to 
model the convection zone. The two isothermal layers are joined through the transition region 
with a steep temperature gradient whose width is ~ 0AH . As a result, the transition region 
physically corresponds to a density contact discontinuous layer. 

For the initial magnetic field we put a magnetic flux sheet in the convectively unstable 
zone and uniform background field. The former represents the flux sheet which produces 
emerging flux and the latter represents the pre-existing magnetic field in the chromosphere and 
the corona. A schematic illustration of the model is shown in Figure 2. 

The magnetic flux sheet in the convectively unstable zone is given as 

~87rp(z) ll/2 



where 



and 



0(z) 



X 



tanh ( 



z ■ 



w 



+ 1 



•tanh 



Zl 



+ 1 



(14) 



(15) 



(16) 



Here, is the plasma beta at the center of the magnetic flux sheet, and z and z\ = z + D 
are the heights of the lower and upper boundaries of the flux sheet, respectively, where D is 
the vertical thickness of the sheet. We use z = —4:H ,D = 2H , /3& = 5 and w = W\ = 0.2H . 
The initial magnetic field strength is ~ 45 (^ 10 3 G) at the center of the sheet. 
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Fig. 3. Vertical distributions of the temperature, magnetic pressure, gas pressure and density in the initial 
condition. The units are 5,600 K for temperature, 6.3 x 10 4 dyn cm -2 for the gas and magnetic pressure, 
and 1.4 x 10 _T g cm -3 for the density, respectively. 

Using the initial plasma beta and temperature distribution mentioned above, we cal- 
culate the initial density and pressure distribution numerically using the equation of static 
pressure balance: 



where B cor is the strength of the background field and 9 cor is the angle of the field from the 
x-axis. We took B cor = 0.3B (75G) and 9 cor = 2n/3. We show the vertical distribution of 
physical parameters in the initial condition in Figure 3. 
A small-density perturbation of the form 



is introduced to the magnetic flux sheet within the finite horizontal domain (— A/4<x<A/4), 
where A is the perturbation wavelength and A is the maximum value of the initial density 
reduction. We took A = 15H and A = 0.1p - The wavelength is nearly the most unstable 
wavelength of the linear Parker instability (Parker, 1966). 

We use periodic boundary conditions in horizontal direction, free boundary conditions 
at the top, and free boundary conditions for B x and B z and symmetric boundary conditions 
for the other parameters at the bottom. 

The size of the simulation box is —70Hq < x < 70H and — 10H < z < 140-f/o? namely 
23,800 km x 25,500 km, which is large enough to include the observed chromospheric jets 
(~ 1,000 km). The grid spacing in horizontal direction is uniform with Ax = O.IHq in — 40i^o < 
x < 4:0H and increases upto 1.5H in x < — 40i^o and 40i^o < x. In vertical direction Az = 0.05H 




(17) 



The uniform background field is given as 
-Bcor = (B cor cos9 cor )x + (B cor sin9 COT )z 



(18) 



Sp = Af(z)p(z) COS (27TX/X) 



(19) 
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in z < 30H and increases up to I.OHq in z > 30H . The total grid number is N x xN y = 828 x 906. 

The numerical scheme we adopted is based on the CIP-MOCCT method (Kudoh et al., 
1999). The scheme can capture the profiles of physical parameters with contact discontinuities. 
Therefore using this method we can investigate the dynamics of the transition region in de- 
tail. We used some of the routines contained in CANS (Coordinated Astronomical Numerical 
Software) maintained by Yokoyama et al. 

3. Numerical Results 

3.1. Evolution of Emerging Flux 

Magnetic loops expand as a result of the convective Parker instability (Parker, 1966; 
Nozawa et al., 1992). The time evolution of the magnetic loops is shown in Figure 4. In 
Figure 4, two jets are indicated by the white arrows (Jet A and B). We focus on these jets and 
investigate how magnetic reconnection creates these jets. 

3.2. Jet Resulting from Magnetic Reconnection near Photosphere: Jet A 

The foot-points of the emerging loops sink into the convection zone because the plasma in 
the lifted flux sheets falls down along the magnetic field and the foot-points become heavy. The 
sinking motion naturally creates an antiparallel field (U-shape structure) and a current sheet. 
This behavior is also found in the models presented by Isobe et al. (2007) and Archontis & Hood 
(2009). As the result of the current sheet formation, magnetic reconnection takes place just 
below the photosphere [z ~ —4Hq) (see Figure 5). In our simulation, the numerical resolution 
is insufficient to resolve the reconnection region or the current sheet. Although we cannot 
capture detailed structures of the reconnection region, we infer that the simulation describes the 
global dynamics like the acceleration of the reconnection outflow. The reconnection rate of the 
numerical reconnection is ~0.1. If a fast reconnection takes place near the actual photosphere, 
the reconnection outflow will generate slow mode waves as found in this simulation. 

The reconnection outflow is drastically slowed down by the dense and high-beta pho- 
tospheric plasma within a short travel distance. The outflow moves upward along a magnetic 
field due to both the inertia effect and squeezing. The upward flow collides with the downward 
flow along the field. As a result, a compressed region is formed, leading to the generation of a 
slow mode wave. This behavior is also found in the model presented by Takeuchi & Shibata 
(2001). 

The slow mode wave shows wave amplification when it propagates upward through the 
gravitationally stratified atmosphere as a result of the law of conservation of energy. Eventually 
the slow mode wave becomes a slow shock. This will be analyzed in detail in Section 4. The slow 
shock produces a chromospheric jet through the interaction with the transition region. Note 
that to guide the slow mode wave generated by magnetic reconnection in the lower atmosphere, 
an ambient field plays an important role. The chromospheric material of the jet is expanded 

10 



Log Density Log Temperature 

10" 8 10" 5 10" 2 10 1 xp 10" 2 3x10" 1 10 1 3x10 2 xTo 




X (170km) X (170km) 

Fig. 4. Time evolution of the density (left) and the temperature (right). The black solid lines show 
the magnetic field lines. The arrows show the velocity vectors, po = 1.4 x 10 -7 g cm -3 , To = 5,600 K, 
C s o = 6.8 km s _1 and r = 25 s. 
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Fig. 5. Jet A: Time evolution of the density (left), the temperature (middle) and the divergence of 
the velocity field (div V; right). The black solid lines show the magnetic field lines. The arrows show 
the velocity vectors. The green solid lines show the transition region. Note that in the images of div 
V the blue and the red regions correspond to the compressed and the expanded regions, respectively. 
po = 1.4 x 10~ 7 g cm" 3 , T = 5,600 K, C s0 = 6.8 km s" 1 and r = 25 s. 
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(see div V distribution of Figure 5) and the density is decreased. This may be important for 
disappearance of jets in images of chromospheric lines. 

We obtained the physical parameters along the field line indicated in Figure 6. Figure 7 
shows the time evolution of the velocity parallel to the magnetic field line, the density and the 
pressure measured along the field line. 

The growth of the shock, i.e., the growth of the amplitude of the velocity parallel to the 
field line, can be seen in Figure 7 (b). The velocity of the jet is — 6C s q ~ 40 km s _1 . This is much 
larger than the local sound speed. Note that the upward velocity of the reconnection outflow 
in the photosphere is 0.6C s o ^ 4 km s _1 . The maximum height of the jet is — 25H ~ 4,300 km 
from the photosphere. The difference between the maximum height and the initial height is 
-17^0^ 2,900 km s" 1 . 

This example shows that even if magnetic reconnection takes place below the photo- 
sphere, tall (— 4 Mm) jets can be produced through the interaction between a slow shock and 
the transition region. It should be noted that the jet is not the reconnection outflow but the 
transition region lifted by the slow shock. 

3.3. Jet Resulting from Magnetic Reconnection in Upper Chromosphere: Jet B 

Magnetic loops push an ambient field and start to reconnect with it as magnetic loops 
expand. Chromospheric material is accelerated after magnetic reconnection. The accelerated 
chromospheric plasma eventually start to rise upward along the magnetic field, becoming a 
chromospheric jet. No hot jet driven by the pressure gradient as in Yokoyama & Shibata 
(1996) coexists with the chromospheric jet. 

Figure 8 shows the time evolution of the density, temperature and divergence of the 
velocity field of the emerging flux region. The density distributions show a shape similar to 
that of observed jets, i.e. a jet with a loop structure (see Figure 1). 

We assumed the anomalous resistivity model in the simulation (the justification of the 
anomalous resistivity model will be discussed in Section 4). As a result, the Petschek-type 
reconnection, in which the standing slow shocks emanate from the magnetic reconnection point, 
occurs between the emerging flux and the pre-existing ambient field. In the following, we call 
the standing slow shocks emanating from the reconnection point the Petschek slow shocks. 

Figure 9 shows the time evolution of the resistive electric current \r]J y \ at the neutral 
point as a measure of the reconnection rate, where the reconnection rate is defined as the 
reconnected magnetic flux per unit time. This is measured in the region — 10H < x < 5H and 
H < z < 10 H . As shown, magnetic reconnection starts at t — 46r. 

Many shocks are formed in the chromosphere. The blue and red regions in div V dis- 
tribution in Figure 8 correspond to the compressed and expanded regions, respectively. The 
shock structures can be discerned as sharp blue regions in Figure 8 (blue regions emanating 
from the reconnection point are the Petschek slow shocks). The enlarged images before and 
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Fig. 6. Jet A: Time evolution of the density (left), the temperature (middle) and the divergence of the 
velocity field (div V; right). The black solid lines show the magnetic field lines. The data shown in Figure 7 
were obtained along the thick black lines on the temperature distributions. Those thick black lines are 
the same field line. 

after one of the Petschek slow shocks crosses the transition region are shown in Figure 10. 

We obtained the physical parameters along the field line shown in Figure 11. The black 
thick field lines are the same field line. Figure 12 shows the evolution of the velocity parallel 
to the magnetic field line, the density and the pressure measured along the field line. The 
enhanced density and the pressure in the outflow region appear as the peaks at distance of 
about 10H in Figure 12 (The positions where the distance is WHq are marked with the black 
points on the field line in Figure 11). We can see that the height of the transition region (i.e. the 
steep gradient in the density profile at 15 H distance from the start point of the line at t = 59r) 
decreases with time during the period between t = 59r and 60r . This is because the reconnection 
inflow reduces the height of the transition region (see Figure 8 and 10). Eventually at t ~ 60r 
the transition region crosses one of the Petschek slow shocks attached to the reconnection point. 

After reconnection started, the inflow and outflow regions are separated by stationary 
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Fig. 7. Jet A: (a) Time evolution of the velocity parallel to the field line, the density and the pressure. 
These data were obtained along the field line shown in Figure 6. (b) A close-up image of the period 
between t = 47r and 51r. SS and TR stand for the slow shock and the transition region, respectively. 
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Fig. 8. Jet B: Time evolution of the density (left), the temperature (middle) and the divergence of 
the velocity field (div V; right). The black solid lines show the magnetic field lines. The arrows show 
the velocity vectors. The green solid lines show the transition region. Note that in the images of div 
V the blue and the red regions correspond to the compressed and the expanded regions, respectively. 
po = 1.4 x 10" 7 g cm -3 , T = 5,600 K, C s0 = 6.8 km s" 1 and r = 25 s. 
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Fig. 9. Jet B: Time evolution of the resistive electric current \rjJ y \ at the neutral point as a measure of 
the reconnection rate. This is measured in the region — lOi^o < x < 5Ho and Ho < z < IOHq. The unit of 
the time is 25 s. 

Petschek slow shocks. We express such a situation in words gthe slow shock is attached to the 
outflow region. h When one of the Petschek slow shocks crosses the transition region, it starts 
to propagate upward. Then the shock is no longer attached to the chromospheric material 
accelerated due to magnetic tension (see the bottom panels of Figure 10) (The propagating 
speed of a slow shock in the low-beta region is almost the same as the local sound speed, which 
is much larger in the corona than that in the chromosphere.) In div V distributions at t = 59. 8r 
and later times in Figure 8, we can see the slow shock (blue region) propagating upward. After 
the interaction, we can find that the steep gradient in the magnetic field strength distribution 
is smoothed out. 

After the interaction between the Petschek slow shock and the transition region, the 
reconnection outflow, which is almost in the horizontal direction, is accelerated, although the 
transition region is not lifted upward yet (see Figure 12 (b)). Considering that the outflow speed 
is comparable to the local Alfven speed, it is worth showing the Alfven speed distribution of 
the reconnection region. Figure 13 shows that the Alfven speed in the outflow region becomes 
large after the interaction. In the region where the Alfven speed increases, the outflow velocity 
also has increased. 

To investigate the interaction process in detail, we performed 1.5D MHD simulations and 
analyzed the interaction between one of Petschek slow shocks and the transition region. We 
summarize the simulation results in Appendix. We found that the outflow is further accelerated 
due to magnetic tension after the interaction between the Petschek slow shock and the transition 
region. When the slow shock crosses the transition region, a slow mode rarefaction wave is 
generated and then propagates into the outflow region. The slow mode rarefaction wave causes 
the density and the magnetic field strength to become low and high, respectively. Thus the 
Alfven speed becomes large (see Figure 13). 



17 



Log Temperature div V 

10" 2 3x1 0" 1 10 1 3x10 2 xTo -1 1 




-15 -10 -5 -15 -10 -5 
X (170km) X (170km) 

Fig. 10. Jet B: Enlarged images of the temperature (left) and the divergence of the velocity field (right) 
distributions before (top) and after (bottom) the slow shock crosses the transition region. The green line 
shows the transition region. To = 5,600 K, C s o = 6.8 km s _1 and r = 25 s. 

When magnetic reconnection takes place, a magnetic field tends to be relaxed and re- 
configured to a relatively straight structure. Reconnection outflows are accelerated as a conse- 
quence of the relaxation process. The Alfven velocity near the outflow region determines the 
outflow speeds. Therefore, the relaxation process will proceed quickly when the Alfven veloc- 
ity increases. Therefore the high Alfven speed leads to quick relaxation of the configuration 
of the magnetic field. The re-configuration process results in the further acceleration of the 
reconnection outflow. 

How is the chromospheric plasma accelerated upward due to magnetic reconnection? 
We noted that the upward motion of the transition region starts around at t = 62r. To clarify 
how the jet is accelerated, we show the forces acting on the chromospheric plasma just below 
the transition region in Figure 14. The plasma in the outflow region feels no strong upward 
acceleration during t = 60r and 62r. The z-component of the velocity shows the gradual 
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Fig. 11. Jet B: Time evolution of the density (left), the temperature (middle) and the divergence of 
the velocity field (div V; right). The black solid lines show the magnetic field lines. The data shown in 
Figure 12 were obtained along the thick black lines on the temperature distributions. Those thick black 
lines are the same field line. The position where the distance of 10 Ho from the start point is marked with 
the black point. 
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Fig. 12. Jet A: (a) Time evolution of the velocity parallel to the field line, the density and the pressure. 
These data were obtained along the field line shown in Figure 11. (b) A close-up image of the period 
between t = 59r and 61r. SS and TR stand for the slow shock and the transition region, respectively. 
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Fig. 13. Jet B: Left: Time evolution of the distribution of the Alfven speed before (t = 58. St) and after 
(t = 60. 2t, t = 60. 8r and 62. Or) the interaction between the slow shock and the transition region. The thick 
white line shows the transition region. The thick black line show the same filed line. The Alfven velocity 
in the region between the dashed lines is increased. C s q = 6.8 km s _1 and r = 25 s. Right: Schematic 
diagrams to describe how the slow shock crosses the transition region. SS and TR stand for the slow shock 
and the transition region, respectively. The cross mark represents the reconnection point. 

acceleration during the period. Note that the chromospheric plasma in the outflow feels no 
strong acceleration due to the pressure gradient force when the outflow collides with the ambient 
field. The chromospheric plasma (reconnection outflow) which is almost horizontally accelerated 
eventually starts to rise upward along the magnetic field after colliding with the ambient field. 

A schematic diagram of the acceleration process is shown in Figure 15. After the outflow 
collides with the ambient field, the outflow plasma moves in the direction along the ambient 
field. Consequently, we obtain the chromospheric jet. This process which we call the "whip-like 
acceleration" is essentially the centrifugal force acceleration. 

The upward velocity (V z ) of the chromospheric jet is ~ 2.5C s o ~ 17 km s _1 . The maximum 
height of the jet is ~ 15Hq ~ 2,600 km from the base of the photosphere (Figure 12). The 
difference between the maximum height and the initial height is ~ 8H Q ~ 1,400 km s _1 . 

We can discern the transition region corrugated at the top of the jet (see Figure 8 and 
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Fig. 14. Jet B: (top) Temperature distribution with the velocity and three force vectors. The velocity 
and the forces are measured just below the transition region along the field line in Figure 11. (middle and 
bottom) Time evolution of the z-components of the velocity and the forces. The unit of the forces per 
unit mass Fq is C s q/t = 0.27 km s -2 , where C s q = 6.8 km s _1 and r = 25 s. 
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Fig. 15. Jet B: Schematic diagram to describes how the the transition region starts to rise upward ("whip- 
like acceleration"). TR and SS stand for the transition region and slow shock, respectively. For a detailed 
discussion on the SS-TR interaction, see Appendix. 
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13). Stone & Edelman (1995) found that parallel slow shocks are unconditionally unstable to 
the corrugation instability. In out simulation, the slow shock crosses the transition region and 
propagates almost parallel to the magnetic field. Therefore the corrugated structure can be a 
result of the corrugation instability. 

4. Discussion 

We studied the acceleration mechanisms of chromospheric jets associated with emerging 
flux using a 2D MHD simulation. It is found that the height of magnetic reconnection deter- 
mines the acceleration mechanism of chromospheric jets. Here we present a schematic diagram 
summarizing the acceleration mechanisms of chromospheric jets. We also give implications for 
observations and discuss the validity of the adopted assumptions. 

4-1- Comparison between One- Dimensional Case and Two-Dimensional Case: Jet A 

We found that a slow mode wave grows to be a shock in the case that magnetic recon- 
nection occurred just below the photosphere. The growth of the shock in one-dimensional cases 
has been studied by many authors (Osterbrock, 1961; Shibata et al., 1982; Heggland et al., 
2007). We investigate how our 2D simulation results differ from the prediction of the idealized 
ID theory in which a flux tube is rigid enough not to be disturbed by the plasma dynamics in 
it and the total energy is exactly conserved in a flux tube (so-called zero-beta approximation). 

Here we review the ID theory. We consider the situation in which a slow mode wave 
propagates upward in a rigid flux tube in a gravitationally stratified atmosphere. In the linear 
regime, we obtain the following relation from the energy conservation of the slow mode wave, 

Apv^C s = constant, (20) 

where A, v\\ and C s are the cross-section of the flux tube, the velocity parallel to the mag- 
netic field and the sound speed, respectively. Considering the magnetic flux conservation 
£M=constant and that C s is almost constant in the chromosphere, we obtain 

^l|£-°- 5 ocp- - 5 , (21) 

where B is the magnetic field strength. In the non-linear regime, this relation is modified due 
to the dissipation processes. Ono et al. (1960) derived the relation in a strong shock case: 

^-^ocp" - 236 (22) 

Note that these relations will hold if the energy is exactly conserved in a flux tube. 

We investigate the dependence of the growth rate v\\B~ - 5 on the density p using our 
simulation results. We show the result in Figure 16 (a) (This figure corresponds to Figure 6 
in Shibata & Suematsu (1982).) In Figure 16 (a), the triangles are obtained from the numer- 
ical simulation. The dependence of the growth rate on the density changes approximately at 
— logp = 2. We separately applied the least square method to the numerical results below and 
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Fig. 16. Jet A: (a) The growth of the slow shock discussed in Section 3.2. Triangles show the numerical 
results. The two solid lines shows the fitting lines calculated using the least square method (one is obtained 
using left 8 points of the numerical data and the other is obtained using the rest points). The dotted and 
dashed lines correspond to the analytic relations (21) and (22), respectively, (b) The distribution of the 
divergence of the velocity field (div V). The black lines are magnetic field lines. The thick black line is the 
line along which the data in (a) were obtained. The light blue lines correspond to the region where the 
plasma beta is unity. The green line corresponds to the transition region. 

above this density value to obtain the fitting lines. These lines correspond to the two solid 
lines. These lines have slopes of 0.36 and 0.15, respectively. The dotted and the dashed lines 
correspond to the analytic relations (21) and (22), respectively. 

Below — logp ~ 2, the dependence of the growth rate on the density slightly deviates 
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from the prediction of the idealized ID theory (21). This is probably because the plasma beta 
is greater than unity in the low chromosphere and the low-beta approximation is broken. This 
means that the energy of the shock wave in a flux tube leaks to other flux tubes. Above 
— logp ~ 2, the discrepancy between the numerically obtained relation and the ID analytic 
relation (22) becomes large. The energy loss of the wave along a flux tube is larger than that 
predicted with the aid of the one-dimensional analysis. This implies that the one-dimensional 
approximation, i.e., the low-beta approximation is broken in the non-linear regime. 

To clarify the cause of the discrepancy, we investigate whether the assumption of the ide- 
alized ID theory are valid. We plot the lines on which the plasma beta is unity in Figure 16 (b). 
It should be noted that the plasma beta behind the shock is close to unity even in the upper 
chromosphere. Therefore the ID approximation (zero-beta approximation) is severely broken in 
the regime and the energy of the slow shock is carried to other flux tubes by fast mode waves. 
If we consider the energy conversion rate from the magnetic energy liberated by magnetic 
reconnection into the kinetic energy of jets, we need pay an attention to this effect. 

4-2. Acceleration of Reconnection Outflow: Jet B 

The reconnection outflow discussed in Section 3.3 is basically accelerated due to the 
whip-like motion of the reconnected magnetic field. The reconnection outflow is accelerated 
after the interaction between one of the Petschek slow shocks and the transition region. Here 
we focus on this point, which has not been discussed in the previous studies. 

We obtain the slow mode rarefaction wave propagating into the outflow region after 
the Petschek slow shock crosses the transition region. Behind the slow mode rarefaction wave, 
the magnetic field strength increases and the density decreases (see Figure 12). Thus the 
Alfven speed increases, as shown in Figure 13. The large Alfven speed leads to the rapid re- 
configuration of the magnetic field. The re-configuration results in the further acceleration of 
the outflow plasma. For a detailed discussion, see Appendix. 

4-3. Schematic Diagram of Acceleration Mechanisms of Jets 

In the case of the jet associated with magnetic reconnection near the photosphere (Jet A), 
the slow mode wave generated as a result of the collision of the reconnection outflow is amplified 
when it propagates upward. The slow mode wave becomes a shock as a result of the law of 
conservation of energy. Through the interaction between the slow shock and the transition 
region, a chromospheric jet is produced. This is essentially the same process as that discussed 
by Shibata et al. (1982), that is, a pure hydrodynamic process. 

In the case of the jet associated with magnetic reconnection in the upper chromosphere 
(Jet B), the chromospheric plasma (reconnection outflow) is accelerated due ti the whip- like 
motion of the reconnected magnetic field. The outflow is further accelerated through the 
interaction between the Petschek slow shock and the transition region. This is an MHD effect 
which cannot be found in the pure hydrodynamic cases (See Section 3.3 and Appendix). 
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Fig. 17. Schematic diagrams of the generation mechanisms of the chromospheric (cool) jet. The spec- 
ulated scenario for the low and middle chromospheric cases (the second scenario from the left) and the 
model found by Yokoyama & Shibata (1996) are also included. Jet A and B are the jets discussed in 
Section 3.2 and 3.3, respectively. 

We self-consistently related the acceleration mechanisms of chromospheric jets with mag- 
netic reconnection under the assumptions adopted in this study. The height of a reconnection 
point is found to determine the acceleration scenario. The acceleration processes are summa- 
rized in Figure 17. We include the model found by Yokoyama & Shibata (1996) to unify the 
acceleration processes and the speculated scenario for the low and middle chromospheric cases. 

4-4- Implication for Observations 

4-4-1 . Implication for recurrent Acceleration of Jets 

Observation found the recurrent acceleration of chromospheric anemone jets (Nishizuka 
et al., 2011). Here we discuss a possible recurrent acceleration scenario in terms of the process 
found in this study. 

We showed that slow mode waves can be generated when reconnection outflows collide 
with an ambient field. Let us assume that the velocity of the reconnection outflows which 
hit the ambient field is time-dependent. This is possible, e.g., if the plasmoids are formed in 
the current sheet, where the plasmoids are the regions of magnetically confined plasma. The 
emergence of the plsasmoids has been theoretically and observationally implied (Singh et al., 
2011; Leake et al, 2012; Singh et al, 2012). 

When the outflow hits the ambient magnetic field, slows mode waves will be generated 
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Fig. 18. Schematic diagram of a recurrent acceleration process of the chromospheric jet as a result of the 
collision of the reconnection outflow. The slow mode waves are generated as a result of squeezing of the 
plasma at the termination region of the reconnection outflow (a mode conversion.) 

by squeezing. This process is essentially the same as a mode conversion process from fast 
mode waves into slow mode waves (e.g. Bogdan et al., 2003). The slow mode waves (or slow 
shocks) will propagate upward and accelerate jets when passing through the transition region. 
Yokoyama & Shibata (1996) discussed the direct acceleration due to the pressure enhancement 
due to the squeezing of the chromospheric plasma by reconnection outflows. Here we note that 
slow mode waves can be generated as a result of the collision of reconnection outflows as in 
Takeuchi & Shibata (2001). This process will be preferable for the recurrent acceleration of 
chromospheric jets. The schematic diagram is shown in Figure 18. 
4-4-2- Relation between Ellerman bombs and Surges 

Observations have suggested that the Ellerman bombs result from the emergence of the 
undulatory flux tubes (Pariat et al., 2004). Watanabe et al. (2011) observationally found the 
morphological evidence for the magnetic reconnection in the deep photosphere. They also found 
that some Ellerman bombs were accompanied by dark Ha surges. Following their study, the 
apparent velocity of the ejections in the photosphere is ~ 10 km s _1 and the Doppler velocity 
of the surge is ^ 40 km s _1 . 

We obtained the tall (~ 4 Mm) jet resulting from the reconnection just below the photo- 
sphere. The reconnection outflow velocity in the photosphere is 0.6C s q ~ 4 km s _1 . The velocity 
of the jet is 6C s o ~ 40 km s _1 , which is much larger than the reconnection outflow velocity. 
These values are consistent with the observations. The shock acceleration scenario can account 
for the difference between the velocity in the photosphere and the velocity of chromospheric 
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jets. 

4-4-3- Ejections from Quiescent Prominences 

We found that the reconnection outflow speed increases through the interaction between 
one of the Petschek slow shocks and the steep density jump layer. Recent observations have 
implied that magnetic reconnection takes place in quiescent prominences, where the quiescent 
prominences are large structure composed of relatively cool (~ 10 4 K) plasma in the corona of 
the quiet-sun regions (Hillier et al., 2011). There are the density jump surfaces between the 
quiescent prominence and the corona. Therefore, if the Petschek-type reconnection occurs in 
quiescent prominences, the reconnection outflow will be accelerated through the interaction. 

4-5. Justification of Anomalous Resistivity Model in Chromosphere 

We adopt an anomalous resistivity model not only in the corona but also in the chromo- 
sphere. The resistivity model is widely used to understand the rapid reconnection processes, 
namely the reconnection rate is close to unity, in a fully ionized and collision-less plasma like 
the corona (e.g. Yokoyama & Shibata, 1996). Using the resistivity model, many observational 
characteristics of explosive phenomena like solar flares and coronal jets have been successfully 
accounted for (for review, see Shibata & Magara, 2011). If the anomalous resistivity model 
is adopted, the resistivity is spatially localized, which leads to the rapid Petschek-type recon- 
nection. Although there is no clear evidence that the Petschek slow shocks are indeed formed 
during reconnection processes, we can utilize the resistivity model to allow fast reconnection to 
occur. 

We should note that the chromosphere is fully collisional and partially ionized. For this 
reason, the characteristics of the resistivity in the chromosphere can be different from that in 
the corona, which makes it difficult for us to predict the reconnection physics there. Therefore 
before Hinode era we had not known whether rapid reconnection processes would ubiquitously 
occur in the chromosphere. Recently Hinode have observed many rapid events including the 
chromospheric anemone jets, which implies that rapid reconnection processes are taking place 
also in the chromosphere. 

It is known that the collision between electron, ion and neutral can play an important 
role in magnetic reconnection (Zweibel & Brandenburg, 1997). Neutrals mainly contribute to 
the gas pressure in a current sheet. The collision between ions and neutrals leads to the diffusion 
of the neutrals from the current sheet. This diffusion, i.e., the ambipolar diffusion, can case the 
current sheet to become thin. When the thickness of the current sheet becomes comparable to a 
microscopic scale, magnetic reconnection will set in. Recently, Isobe (in prep) showed that even 
with uniform resistivity, the localized (or non-uniform) neutral distribution in the current sheet 
induces the local thinning of the current sheet, leading to the Petschek-like reconnection. This 
scenario may account for fast reconnection processes in the chromosphere. Because there is a 
possibility that the Petschek-type reconnection takes place in the chromosphere, we adopted the 
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anomalous resistivity model to obtain it. We believe that we can mimic reconnection processes 
in the partially ionized and fully collisional plasma using this resistivity model. 

4-6. Effect of Heat Conduction 

Generally speaking, the heat conduction weaken shocks because the heat energy liberated 
at the shock fronts can escape in front of the shock surfaces. We do not include the effect of the 
heat conduction in this paper. We consider whether or not the heat conduction significantly 
affects the simulation results. 

The timescale of the heat conduction can be estimated as 

Tcand ^ 3 X l^n^T^' 2 S, (23) 

where n e , /, X are the electron density, typical length scale, and the temperature, respectively. 
The heat conduction is important at a shock front because the temperature gradient is steep. 
In our simulation, the minimum length scale is the smallest grid size. By taking I = 8.5 km (grid 
size), n e = 10 12 cm -3 and T = 10 4 K, we estimate the heat conduction time r con d at <~ 2 x 10 4 s. 
Considering that the propagation time of shocks in the chromosphere is less than 100 s in our 
simulation, we can neglect the effect of the heat conduction on the shock propagation in the 
chromosphere. 

Heating due to the thermal conduction and cooling due to the radiative loss control the 
location of the transition region. In slow shock acceleration scenario, the maximum height of a 
chromospheric jet is a function of the height of the transition region before a slow shock comes 
(Shibata et al., 1982). Therefore the maximum height of jets will be affected if we include the 
thermal conduction. However, the acceleration scenario in which slow shocks lift the transition 
region will not be changed. 

The heating and cooling processes will affect the coronal EUV emission. Therefore one 
should include the heat conduction as well as cooling terms into the energy equation if one tries 
to directly compare results of numerical simulations with observations like De Pontieu et al. 
(2011) (see e.g. Heggland et al., 2009). 

We thank Dr. A. Hillier, Dr. S. Imada and Dr. Y. Katsukawa for fruitful discussion. 
We also thank Dr. E. Priest for useful comment. Hinode is a Japanese mission developed 
and launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC (UK) 
as international partners. It is operated by these agencies in co-operation with ESA and NSC 
(Norway). This work was partially supported by the JSPS Core-to-Core Program 22001. This 
work was supported by the Grand-in-Aid for the Global COE Program "The Next Generation 
of Physics, Spun from Universality & Emergence" from the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT) of Japan. HI is supported by the Grant-in-Aid for 
Young Scientists (B, 22740121). The numerical calculations were carried out on SR16000 at 
YITP in Kyoto University. The authors thank the Yukawa Institute for Theoretical Physics at 



30 



(a) Corona 



Chromosphere 



Transition 
region 




(k) Reconnected magnetic field 

^C^^^^ Density 

"''^^^^ jump layer 




Reconnection 
point 



Fig. 19. Jet B: (a) the 2D situation discussed in Section 3.3. (b) Simplification of the problem in which 
the density jump layer is located in the inflow region. 

Kyoto University. Discussions during the YITP workshop YITP-W-12-08 on "Summer School 
on Astronomy & Astrophysics 2012" were useful to complete this work. 

Appendix 1. 1.5 Dimensional Study on Interaction between Slow Shock and 
Transition Region 

The process found in Section 3.3 can be considered as a magnetic reconnection problem 
with a density contact discontinuous layer (transition region) in the inflow region (see Figure 19). 
We investigated the interaction process between the Petschek slow shock and the transition 
region. 

A. 1.1. Basic Equations and Assumptions 

We simulated the time evolution of the reconnected magnetic field using 1.5D MHD 
simulations. The gravity is neglected, therefore the results here are scale-free. The basic 
equations are as follows: 
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Fig. 20. The initial condition of the 1.5D simulation in the case that the ratio of the density in the 
chromosphere (high-density region) to the density in the corona (low-density region) is 100 and the plasma 
beta is 0.1. 
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Table 1. 


Parameters 




Parameter 


Value 




Pch/pcr 


10, 100, 1000 


P 


0.01, 0.1, 


1, 10 


[degrees] 


15, 20, 25 


, 30 



x p i y td 

c y c 



(A6) 
(A7) 



A. 1.2. Initial and Boundary Conditions 

Figure 20 shows the initial condition of the simulation in the case that Pch/pcr is 
100, where pcu and pcr are the density in the chromosphere and the density in the corona, 
respectively. The initial angle of the reconnected field line 0, the plasma beta (3 and the ratio of 
the density in the chromosphere and the density in the corona Pch/pcr are parameters. The 
parameters used are summarized in Table 1. We investigated 48 cases. 
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The initial distributions of the physical parameters are as follows: 
P = Pch ~ (pch ~ Pcr)\ (tanh (- — — ) + l) 




(A8) 
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where B y Snpo/fi, w 0.02, x c = 4, x ( i = x c + 0.6, pcTf 1 and p = I/7. The chromosphere 
and the corona correspond to the regions of < x < 4 and 4 < x, respectively. The velocity is 
normalized by the sound speed in the chromosphere. 

The size of the simulation box is < x < 240. The grid spacing is uniform with Ax = 
0.0027 in < x < 10 and gradually increases in 10 < x. The total grid number is 8200. We use 
free boundary conditions. The scheme used is a CIP-MOCCT scheme. 

A. 1.3. Numerical Results 

In Figure 21, we show the time evolution of physical parameters for the case that 
Pch/pcr = 100, /5 = 0.1 and 9 = 20 degrees. We can discern two slow shocks are attached 
to the outflow region until t ~ 0.6. They corresponds to the Petchek slow shocks. The transi- 
tion region (TR) is advected by the inflow, which corresponds to a fast mode rarefaction wave, 
into the outflow region. After one of the slow shocks crosses the transition region, the slow 
shock (SS) propagates in the positive x-direction. Simultaneously, a slow mode rarefaction wave 
(SR), which corresponds to the reflection wave, is generated and propagates in the negative 
x-direction. 

From Figure 21, we find that the outflow velocity increases after the interaction between 
the slow shock and the transition region. We obtained a relation between the outflow velocity 
before the interaction ^outflow 1 and that after the interaction ^ utflow2: ^outflow2 ~ 2^ ou tflowi (see 
Figure 22). 

A. 1.4- Interpretation of Numerical Results 

Why does the outflow velocity increase through the interaction between the Petschek 
slow shock and the transition region? We note that the slow mode rarefaction wave propagates 
into the outflow region after the interaction. Behind the slow mode rarefaction wave, the 
magnetic field strength increases and the density decreases. Thus the Alfven speed increases. 
The increase of the Alfven speed can be found in our 2D simulation (see Figure 13). This leads 
to the rapid relaxation of the configuration of the magnetic field. The re-configuration of the 
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3 4 53 4 53 4 5 
X X X 

Fig. 21. Time evolution of the velocity in x-direction, the velocity in y-direction, the gas pressure, the 
density, the magnetic field strength and the Alfven speed. TR, SS and SR denote the slow shock and the 
slow mode rarefaction wave, respectively. 
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Fig. 22. Relation between the outflow velocity divided by the x-component of the unperturbed slow mode 
speed in the chromosphere before the interaction and that after the interaction, pch and pcR are the 
density in the chromosphere and the density in the corona, respectively. The diamonds, triangles, squares 
and crosses are the results of the cases with /3=0.01, 0.1, 1 and 10, respectively. The four points of each 
mark are the results of the cases with #=15, 20, 25 and 30 degrees, respectively. 
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Fig. 23. Time evolution of the reconnected magnetic field. TR, SS and SR denote the transition region 
(density contact discontinuous layer), the slow shock and the slow mode rarefaction wave, respectively. 

magnetic field causes the outflow plasma to be further accelerated. The schematic diagram is 
shown in Figure 23. 

Considering that the outflow plasma is further accelerated due to magnetic tension after 
the interaction, we estimate the outflow velocity after the interaction as follows: 

<%~h B -VB (M5) 

A^J^L^At (A16) 
4np ox 

4*pWv su ,' [ ' 

where Av is the increase of the outflow velocity through the interaction, W is the length scale 
(see Figure 23) and v$m x is the slow mode speed in the x-direction in the outflow region before 
the interaction. In the outflow region where the magnetic field is in the x-direction, 
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vsu*={ C ' ( f <2/ :\ (A18) 
V Ax W>2/ 1 ), 



where V Ax = B x /y/4wp is the Alfven speed in x-direction. For all the cases, the plasma beta in 
the outflow region is greater than 2/7. Thus vsmx = V Ax . Therefore 

a B x By W 

Av~- — -f- (A19) 

4np W v S mx 

B x B V W 

(A20) 



^pWV Ax 

~7=F= ( A21 ) 
V Ay , (A22) 



where V Ay = By/^Anp is the Alfven speed in y-direction. Before the interaction, the reconnec- 
tion outflow velocity iWflowi is V Ay . After the interaction, -u ou tflow2 = V Ay + At; ~ 2 V Ay — 2-u outfl0 wi • 
This relation shows no strong dependence on the plasma beta and the ratio of the density jump, 
which is consistent with the numerical results. 
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